Error correction in short time steps during the apphcation of quantum gates 
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We propose a modification of tlie standard quantum error-correction metiiod to enable tlie cor- 
rection of errors that occur due to the interaction with a noisy environment during quantum gates 
without modifying the codification used for memory qubits. Using a perturbation treatment of the 
noise that allows us to separate it from the ideal evolution of the quantum gate, we demonstrate 
that in certain cases it is necessary to divide the logical operation in short time steps intercalated 
by correction procedures. A prescription of how these gates can be constructed is provided, as well 
as a proof that, even for the cases when the division of the quantum gate in short time steps is not 
necessary, this method may be advantageous for reducing the total duration of the computation. 
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I. INTRODUCTION 

Since the beginning of quantum-computation theory, 
it has been known that decoherence and other forms of 
external interference in the quantum bits (qubits) pose 
a difficulty in implementing real quantum computers [ij . 
Various methods have been devised to cope with these 
errors, including the quantum error-correction theory, 
which consists of encoding quantum states of k logical 
qubits in n > k physical qubits, thus generating a redun- 
dancy of information that can be used to identify and 
correct errors Q . Many of these codes have since been 
designed for different kinds of quantum channels. 

It has been shown that quantum codes can be em- 
ployed in quantum circuits as long as the gates arc fault- 
tolerant The traditional fault-tolerant methods con- 
sider that the errors can be described as happening after 
the physical gates are applied [3| , which is a harmless as- 
sumption when dealing with a code capable of correcting 
any generic error up to a certain number of qubits. In 
this case, any kind of error that may emerge from the si- 
multaneous interaction of the gate and the environment 
will still be correctable. 

This assumption, nevertheless, does not remain true if 
we take into account quantum gates that have finite-time 
durations which are not negligibly small in comparison 
with the decoherence time Q. When we let both the 
errors and the gates occur simultaneously, unexpected 
new kinds of errors may affect the encoded qubits, which 
may render the correction procedure impractical if we 
are using a code that was designed to correct the noise 
of a specific quantum channel. Given that such channel- 
adapted codes can be more efficient than generic ones @ , 
a method of correcting errors that occur during the gate 
using the same specific code developed for the case of a 
memory qubit deserves attention. 

In this article, we propose an alternative method of 
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correcting errors caused by an external environment oc- 
curring during a logical operation that consists of repeat- 
edly applying the quantum error-correction procedure 
for memory qubits. Employing the analytic theoretical 
methods of standard quantum coding theory and of the 
dynamical evolution of a statistical quantum ensemble 
given by the von-Neumann equation, we lay down the 
conditions for the repeated error correction to be neces- 
sary. We also show how it can be applied to correct a 
universal set of quantum gates, exemplified in the case of 
the threc-qubit phasc-crror-corrccting code |7[ • 

Special attention is given to the form how such gates 
must be constructed. While fault-tolerant circuits em- 
ploy transversal gates, which have been proven incapable 
of constructing a universal set for codes that correct a 
general one-qubit error Q , our method bypasses such re- 
strictions by employing gates which cannot be factored 
in a tensor product of individual-qubit operations, being 
therefore non-transversal. The use of such gates does not 
cause additional errors with significant probability, but 
may require some special accomodations, as explained in 
the body of the article. 

This article is structured as follows: in Sec. II, we ex- 
amine the essential mathematics necessary to understand 
the small-step method. In Sec. II, we explain how the 
method works, and in Sec. IV we show how the gates can 
be implemented within this framework. Sec. V is dedi- 
cated to describe how the gates can be implemented by 
an approximative method when multiplc-qubit interac- 
tions are prescribed. We conclude in Sec. VI, presenting 
some further perspectives of work. 



II. MATHEMATICAL FORMULATION 

In the description of a quantum computer prone to 
errors, we consider that the qubits are part of a sys- 
tem S that interacts with an environment E, usually 
modeled as bosonic baths. Even though we are inter- 
ested only in the reduced density matrix of the qubits 
ps{t) ~ T^TE {psE{t)}, we must take into account the 
complete system, described by the total density operator 
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PS sit), if we want a unitary evolution given by the von 
Neumann equation: 



■^PSEit) = -i [H,psE{t)] , 

where H is the Hamiltonian operating on both system 
and environment and we chose natural energy units, so 
that h ^ 1. 

The solution of this equation is found to be 



PSE 



(1) 



where the time-evolution operator UsE{t) is e~*^* if the 
Hamiltonian is time- independent. This Hamiltonian can 
always be split in three terms: two that act on the system 
and environment separately {Hs and He) and one term 
of interaction that acts on both (i?int)- The latter must 
be weak enough for the errors to be rare - otherwise, 
the recovery of the original state will be impossible. To 
make this assumption explicit, we write the interaction 
term multiplied by a small number A: 



H = Hs 



To treat the evolution due to the interaction with the 
environment as a perturbation in respect to the case of an 
error-free system, we expand the time-evolution operator 
UsE{t) in a power series of A: 



UsE{.t) = Us{t)UE{t) 



l + X]A"C„(t) 



(2) 



where Us{t) and UE{t) represent the evolution due to the 
system and environment Hamiltonians alone. The C„(t) 
are operators that can be found by: 



Cnit) = uimlit) 



1 a" 



■UsE{t) 



(3) 



A=0 



Replacing the expansion of the operator UsE{t) given 
in Eq. ^ in the final state of the system from Eq. ([IJ , 
and taking the trace over the environmental degrees of 
freedom, we find, up to the second order in A: 



Clit) = -Ci{t), 

Applying Snider's identity [§| to calculate the deriva- 
tive of the time-evolution operator, we find: 



g^UsEit) - ~lUsE{t) 



dt' UlE{t')H,^tUsE{t') 



This equation allows us to calculate both Ci (t) and C2 (t) 
through Eq. ©: 



Ci{t) = -i I dt' ul{t')Hi{i')Us{t'), 



C2{t) = Mt)C,{t)~- / dt' 



dt" 



ul{t')Hj{t')Us{t'),Ul{t")Hi{t")Us{t") 



where Hi{t) = U^^{t)HintU E{t) is the interaction-picture 
Hamiltonian for the memory, when there is no gate acting 
on the qubits. 

Therefore, the final state of the system will be: 



ps{t) = Us{t) I dt'TYE[-i\[Hi{t'),psE{0) 



--A^ / dt' 



[Hi{t')Hi{t"lpsE(S))] 



dt" 



Hi{t'),Hi{t") ,pse{0) 



+\'J^ dt"Hj{t')psE{0)Hi{t")j Ulit) 

+Us{t)psi0)ul{t) + O{X^), 

where we defined Hi{t) = ul{t)Hi{t)Us{t) as the 
interaction-picture Hamiltonian for the case when the 
gate is being applied. The first-order term is usually 
neglected [lO[, which can be achieved by the imposition 
that Tr^; {if/(t')p£;(0)} = 0. This condition is especially 
valid when the environment is a bosonic bath, in which 
case the interaction Hamiltonian usually contains terms 
with a single annihilation (0^) or creation (a^) operator: 



ps{t) = us{t)psmut) 

+XUs{t)TTE{[Ci{t),psEm} ulit) 
+X^Us{t)Ti-E {[psE{0)Ci{t),Cim ulit) 
+X^Us{t)TYE {[Chit), PsEiO)]} ulit) + 0(A3), 

where we imposed the condition that UsEit) is unitary, 
so that: 



Hint = 2_^Skak -l-h.c, 

k 

where the 5*^ represent operators that act on the system. 
As the trace of one creation or annihilation operator al- 
ways vanishes: 

Tr£{afc} = Tr£{4} = 0, 



3 



the partial trace over the whole interaction Hamiltonian 
will be zero. 

Considering initially the state of the system separable 
from the environment - before t = no errors were oc- 
curring, which means that S and E cannot be entangled 
- we find the following equation for the final state: 



psit) = ^Xh^Us{t)TvE{£s{t)psEmulit) 
+Us{t)ps{0)uUt)+O{X^), 



(4) 



where the error superoperator £s{t) acts on the density 
matrix according to: 



Hiiut), p,Hi{vt) 



dv 



Hi{ut),Hi{vt) 



The second term on the right-hand side of Eq. Q rep- 
resents how the system evolves when it is subject only to 
the gate Hamiltonian Hs- Therefore, the superoperator 
£s{t) represents the errors that occur, with a probabil- 
ity proportional to X^t'^, while the gate is being applied. 
However, this superoperator depends on the type of the 
gate. An arbitrary quantum code should not be designed 
to fix every kind of error, just those that occur with a 
probability of the order of when the qubit is in the 
memory. As our objective is to avoid the impractical 
case where a different code is used for each gate applied, 
the errors that occur during the gate must be of the same 
type as those that occur. The latter will be denoted by 
the superoperator £o{t), which is the same as £s{t), but 
with the gate Hamiltonian Hs set to zero: 



£o{t)p = £s{t)p\Hs=o 



du\ / dv[Hi{ut),[p,Hi{vt 



dv[[Hi{ut),Hi{vt)],p] 



III. THE SHORT TIME STEP METHOD 

A sufficient condition for the errors during the quan- 
tum gate to be correctable is the commutativity of the 
gate operator Us{t) and the interaction-picture Hamil- 
tonian Hi(t). In this case, the interaction Hamiltonian 
is identical to the case without gate {Hj{t) = Hi{t)), 
so that £s{t) ~ £o(t) and the errors that occur are of 
the same type that happen in the memory. This fact is 
also evident if we expand the time-evolution operator in a 
part that represents the gate and another that represents 
the errors: 



UsE{t) 



-i(Hs+HE + Hi„t)t 



Us{t)e 



-i(HE+Hi„t)t 



if [Hs,Hiat] = 0. In this case, the evolution of the qubits 
due to the gate is separable from the evolution due to the 
noise. Therefore, there is no need to modify the standard 
error-correction procedure. 

However, when those two terms do not commute, the 
kinds of errors generated by £s{t) may not be correctable 
by the code. To solve this problem without having to re- 
sort to more resource-consuming codes, we propose the 
application of the quantum-gate operation of total du- 



ration tg in N short time steps At = N ^tg. After 



each of these steps, the standard error-correction proce- 
dure is applied, until the system reaches the desired final 
state. The ordinary error-correction procedure works in 
this method because, for small time intervals, the Hamil- 
tonians Hs and Hiit) commute approximately, resulting 
in an interaction-picture Hamiltonian identical up to the 
zeroth order in N to the case when no gate is being ap- 
plied: 



Hi{At) = Hi{At) + O(iV-i). 

As each term of the error superoperator contains two 
expressions of the kind Hj{At), the final result for the 
expansion of £s{At) is: 



£s{At)p ^ dul^J^ dv[Hj{uAt),[p,Hj{vAt)]] 
dv [[Hi {uAt) , Hi {vAt)] , p 



+0{N-^) 
= £o{At)p + 0{N-'). 

which leads to the following state of the system after the 
short time interval, according to Eq. Q: 



psit + At) 



2j-2 



XH 

2A^2 



Us{At)TrE {£o{At)psE{t)} ul{At) 



~UsiAt)ps{t)ul{At) + 0(A27V-3). (5) 



Choosing the number of steps N as approximately 
the probability of finding an uncorrectable error 
is of the order of X'^N~^ ^ A^/^, negligible in comparison 
with the probability of finding a correctable error, which 
is represented by the first term in the right-hand side of 
Eq. (O, which is of the order of A^iV"^ ~ A^. Therefore, 
the correction method designed for memory qubits can 
be applied during the gate with low probability of failure. 

The repeated application of the error correction does 
not increase the probability of an uncorrectable error, 
originally of the order of A'^: if the gates are split in 

parts, each with a probability X^N~^ of suffering an 
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uncorrectable error, the final probability of error is of the 
order N (A^7V~^) ~ A'^, the same order as in the memory. 

This method also does not increase the duration of the 
quantum-gate operation: making the reasonable assump- 
tions of instantaneous measurement and instantaneous 
processing of the syndrome, the probability of having to 
apply a correction procedure will be of the order of if 
the gate is not divided in smaller steps. When we apply 
the correction procedure in N small time steps, the prob- 
ability of having to apply a correction is of the order of 
N (AiV-i) ~ A^/2^ Therefore, the new average duration 
of the gate, which is dependent on the number of times a 
correction procedure has to be applied, is actually shorter 
than the standard method. 

The only difficulty we have not dealt with so far are 
the errors that may occur while the correction gate it- 
self is being applied. However, one must note that the 
probability of a complete gate requiring a correction pro- 
cedure (estimated above as of the order of A^7V~^), and 
the correction procedure itself suffering some sort of er- 
ror (A^, as in the case of an ordinary unprotected gate) 
yields a result of the order of A^iV~^ ~ A®/^. This is 
much less than the probability of an uncorrectable error 
(given above in A'^), and, therefore, safely negligible. 

It is important to note that even in the case when the 
partition of the gate in small steps may not be necessary 
(as when the logical gates all commute with the interac- 
tion Hamiltonian or when the unexpected errors are for- 
tuitously correctable by the code), this method can still 
be useful to reduce the probability of an error happening 
and, therefore, the overall duration of the computation. 
Indeed, as demonstrated above, a quantum-gate opera- 
tion with a probability of approximately A^ of having to 
be corrected, with this method is only required to be 
corrected with a probability of the order of A^/^. As the 
necessity of error corrections is smaller, the computation 
may be performed faster, and the errors that occur while 
the correction gate is being applied can be safely ignored. 

Additionally, we may note that this method does not 
require the excessive partition of the gate. If the error 
is of the order of A^ ^ 10^"*, then A ~ 10^^, and our 
prescription requires only iV ^ 10 time steps. However, 
it is possible to increase the number of time steps to 
further reduce the probability of any error occurring. The 
decrease in the probability of error is only limited by 
experimental capability of partitioning the gates. 

The caveat, however, is that our proposed method re- 
quires a gate Us{t) that keeps the state of the qubits 
inside the quantum code for every instant t, not only 
for the final instant tg. Otherwise, the syndrome mea- 
surement will detect errors when none have happened. 
However, the construction of logical gates that respect 
this restriction does not pose an unsurmountable obsta- 
cle if we employ non-transversal gates explained in the 
next section. 



IV. GATES 

Logical gates are intended to reproduce, on an encoded 
logical state, the effects of some physical gate G on a 
physical qubit state Identifying the logical gates 

and states by the subscript L, if the physical gate acts 
according to G\iIj) = then the equivalent logical 

gate must act according to G l ~ W) l- 

All n-qubit quantum gates can be reproduced by some 
combination of CNOTs and single-qubit rotations fllj . 
The correct construction of these two gates, therefore, is 
sufficient to render any computation protectable by this 
method. Any rotation gate can be characterized by three 
angles 0, 0, and 

Ui-ot{(^, 4'^ v) — cos ip — ism (p 

X [az cos 9 + sin 9 cos (j) + ay sin 9 sin </>] . 

The CNOT gate, on the other hand, has always the same 
structure: 

C/CNOT = (^^f^) + (^^) ^2.., 

where the index 1 represents the control qubit, while 2 
represents the target. 

The Hamiltonians used to generate these two gates are, 
for the rotation: 

Hiot{(^i ^) = '^o i'^z cos 9 + sin cos — icTxCz sin sin ^) , 

(6) 

which must be applied for a period tg = cjq ip, and, for 
the CNOT, 

i?CNOT = Wo (—^^^ ' ^"^^ 

where the end time must be chosen so that tg = loq^t:. 

To create the equivalent logical gates, we must replace 
the physical bit flips and phase flips cr^ by their logical 
equivalents in Eqs. (jG]) and ([2]). In this case, the corre- 
sponding time evolution of the logical gate will retain 
the form of the physical gates as long as the logical Pauli 
matrices are Hermitian and still satisfy the identities 

{aL,.f = {<JL..f = 1, (8) 

'7L,z<7L,x = —(^L,x(^L,z- (9) 

The time-evolution operators for instants of time before 
the gate operation is complete will be combinations of 
the logical <Jl,x, ol,z and identity operators. As aL,x 
and aL,z are constructed with the intention of keeping 
the states inside the code subspace, both the rotation 
and the CNOT are guaranteed to keep the qubits inside 
the code for any instant of time. 
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Therefore, we just need to know how to define the log- 
ical Pauli matrices that satisfy Eqs. ^ and We il- 
lustrate the construction of them for the triple-repetition 
phase-error-correcting code, defined as: 

HL = (8)f45E(-ir"NV (10) 

k=l v V ^ „=o / 

For this code, it can be easily seen that (Ji ^ fulfills the 
role of logical phase fiip, and o-i,z<^2.zO'3.z, the role of bit 
flip. As these operators belong to the Pauh group, it is 
clear that they satisfy Eqs. ^ and as required for 
the method to be applicable. In this case, the partition 
of the gate in short time steps is also required, because 
the logical phase-flip gate ctl^z does not commute with 
the interaction-picture Hamiltonian for the phase error 
channel, which contains the Pauli matrix az- 

A final concern that must be addressed is that this 
method of constructing gates precludes the correction of 
errors due to an incorrect duration of the gate, besides 
requiring the use of Hamiltonians that may involve the 
interaction of three qubits. These problems can be cir- 
cumvented if the gate itself is also simulated through the 
short-time application of simpler Hamiltonians, as ex- 
plained in next section. 



so that it is expected that a sequence of such operators 
will be capable of effectively simulating the unphysical 
gate we wish to apply. 

As these two matrices do not commute, we can expect 
that application in sequence of e~*'^i-'^2'="'^e~"^^''''^^-^^ and 
the same couple of operators with the sign of e inverted 
will not yield a result proportional to the identity. In- 
deed, up to the second order in e, what we actually have 
is: 



g — io'l.j(T2,xeg-iCT2,afT3,jegi(Ti^j(T2,x£gio-2,j,f3,je 

l-2ie^ai,z<72,z'J3,z + 0{e^), (11) 

which is equivalent to an approximation of the three- 
qubit Hamiltonian up to the second order in e. A numer- 
ical comparison of this simulation with the ideal non- 
physical gate is presented in Fig. [TJ where it can be seen 
that the order of the magnitude of the error will be 10"'^ 
if the number of steps N approaches 10^. 

To improve on the quality of the simulation, we can 
expand Eq. pip up to the third order in e, finding: 

g-iO'l,,0-2,iCg-io-2,H''"3,iegiO'i,jO-2,a;CgiCT2.y(T3.-C _ 

1 - 2ie^ai,z(T2,zCr3,z + 2ie^ (cr2,y(T3,2 - ai^z(^2,x) + 0(e''). 



V. SIMULATING THE GATES 

It can be seen that some of the gates needed by the 
method described in previous section include three-body 
interactions, which require non-physical Hamiltonians. 
To counter this problem, we can divide each of the 
multiple-body interactions in a series of two-body gates, 
which is valid for small time steps: 

g»e<Ti,3iO-2,32---<^->,j„ =Y\_e"'''''"k-'k'^"k-'k +0{e'^). 

k 

The method used here to simulate the multiple-body 
Hamiltonians, inspired by those already developed for 
long periods [12| , will take advantage of the fact that our 
gates arc applied to small time-steps, so that the approx- 
imation given above is reasonable. 

This scheme will be illustrated for the case of the three- 
qubit code from Eq. ([TU)) . in which case the aL,x gate re- 
quires the simulation of a three-qubit Hamiltonian, while 
the CNOT requires a four-qubit one. In the first case, 
the time-evolution required is of the form e~'^"''^-'"''^-'"'''-''^, 
which must be split into a product of exponentials that 
contain only two Pauli matrices. For this purpose, we 
wiU use the operators e~*'^i'='^2'="'^ and e"*'^^"'^^-^^, cho- 
sen because the product of their arguments retrieves the 
Hamiltonian we are trying to simulate: 

{f^l.,zCr2,x) (CT2,yO'3,z) — i<^l,z(^2,z'y3,z, 



The additional term on the right-hand side can be turned 
into two physical exponentials (containing only two Pauli 
matrices) and put together with the other two-qubit 
terms, so that the expansion up the third order becomes: 

g-2i£^(Ti,^CT2,iCT3,i _ g — i<'"i^j;Cr2,xeg-io'2,aCr3^je 
y g*fl,i<^2,xegio'2,HCr3 

g - 2ie^ 0-2 , H , i g2i<:^ <Ti , J <T2 , X 

+0(6^). (12) 

The precision of this new simulation is also portrayed 
in Fig.[Tl showing that the threshold of 10"'' probability 
of error is achieved with just N ^ 10^ steps. 

The logic gate CNOT requires a four-qubit operator of 
the form e~"^i'^'^2,„cr3,jiT4,xe^ ^j^g gj-g^ expansion of which 
can be simply obtained by adapting Eq. ([T^ : 

g — 2ie^cri^jCr2,iO"3,iO"4,rc _ ^ — 101,^02, ^as^^t ^-i(73^y(y4,_^t 

j^g*fl,i<^2,i<T3,xegio'3,„(T4^xe 

-^g-2ie^cr3 „cr3,xg2ie^cri^^cr2,jtT3^x 

+0(6^). (13) 

The problem now is to simulate the matrices of the 
type e~"^i'^''2''°'3'^^, which are also non-physical. In order 
to keep their errors within the limits of this expansion, 
which correspond to the fourth order in e^, an expansion 
up to the seventh order is required: 
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FIG. 1: (color online) Fidelity of the logical gate for the 
three-qubit code, after split in physical two-body Hamiltoni- 
ans. The two upper curves represent the two less efficient 
second-order approximations, while the two lower represent 
more precise third-order results. Both approximations were 
tested for two initial states, which correspond to the computa- 
tional basis of the code subspace. The y-axis is in logarithmic 
scale to better represent the fast increase of the fidelity. 



FIG. 2: (color online) Fidelity of the logical gate CNOT for 
the three-qubit code, after split in physical two-body Hamil- 
tonians. The different curves represent four different initial 
conditions, which are the computational basis for two encoded 
qubits. The difference between the two curves is due to the 
fact that, depending of the value of the control qubit (sec- 
ond qubit), the operation on the first qubit is distinct. Both 
curves have the same order of magnitude, nevertheless, as can 
be seen from the logarithm scale. 



„i(2£=--8c''/3-56eV45) 



-ieCTi,-cr2,:c i£cr2,ycr3,^ i(6£^ + 16e''/5)(Ti.j(T2.i 



^g-i(2£S-56eV5)<T2,y(T3.:.^ 
Xg-2ie3<T2,ya3.. +0(e8)^ 



-2ie (7i^z<72,'. 



Using the results of Eq. ([H)) in conjunction with those 
of Eq. p3p . we obtain a simulation of the CNOT that 
reaches the threshold of 10^'* for the probability of er- 
ror when the number of steps is of the order of 10^, as 
illustrated in Fig. [51 

Should the experimental system require, the precision 
of the gates can be improved even more by the expan- 
sion of the simulation to even higher orders of magnitude. 
Likewise, it is not straightforward to repeat this proce- 
dure for a higher number of qubits. In so doing, we are 
not only able to simulate the multiple-body gate with 
two-body interactions, but we employ a series of opera- 
tors that do not preserve the code subspace, rendering 
errors in the duration of the gate again detectable. 

An additional aspect that can be considered is that 
this method of simulation need not be necessarily applied 
only in this context. When quantum simulations may be 
necessary, this method could be applied independently 
of the protection scheme to reproduce the effects of non- 
physical Hamiltonians. 



procedures can protect against external errors that hap- 
pen during the quantum gate, even for a code designed for 
a specific quantum channel. This is achieved by decreas- 
ing the probability of errors when these are incompatible 
with the code under normal circumstances. The only re- 
quirement of this method is the possibility of defining 
a set of logical Pauli matrices for the code that satisfy 
Eqs. ([8|) and (|9]). We have shown this to be possible for 
the simple three-qubit phase-error-correcting code, but it 
is expected that it can be adapted to any quantum code, 
as Eqs. ([S]) and ^ are derived from the operational def- 
initions of the Pauli matrices. 

Moreover, in this process we have demonstrated how 
the small-step method can be used to reduce the proba- 
bility of even correctable errors, thus increasing the speed 
of the computation; and have presented a simulation 
method for the non-physical Hamiltonians that can be 
used independently of the error correction. 

The main phenomenon of protection of quantum codes 
by repeated measurement, however, is not entirely un- 
expected: a syndrome measurement performed in close 
intervals restricts the probability of an error happen- 
ing [l^. Research by our group l3| suggests that the 
limit of measurements performed continuously while the 
errors occur may protect even more effectively the qubits. 



VI. CONCLUSION 

In conclusion, we have demonstrated that applying the 
quantum gate in small steps intercalated by correction 
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